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Abstract 
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arbitrary "^^^^Lj angular-spin states are determined at next-to- leading order in the nonrelativistic 
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I. INTRODUCTION 



The use of techniques and concepts from effective field theories (EFT's) has lead to an 
encouraging record of achievements for the description of nonrelativistic fermion-antifermion 
systems. The formulation of a nonrelativistic EFT for such systems has gone through a 
number of important conceptual states Bodwin, Braaten and Lepage [2I, 0] initially 
formulated a method for separating fluctuations at short distances of order the heavy particle 
mass m from those at long distances related to the nonrelativistic momentum and energy, 
mv and mv"^ {v being the relative velocity) and the hadronic scale Aqcd- Subsequent work 
helped to clarify the power counting in v, and the relevant degrees of freedom needed to 
construct an EFT. 

For very heavy quarkonium systems, characterized by the scale hierarchy m ^ mv S> 
mf ^ > Aqcd the formulation has reached a mature state. In Ref. [1] it was realized that it is 
necessary to distinguish soft mv) and ultrasoft (~ mt>^) fluctuations and in Ref. [1] that 
the original NRQCD action [3| is problematic concerning the simultaneous power counting 
of soft and ultrasoft terms. In particular, the Lagrangian needs to be multipole-expanded 
in the presence of ultrasoft gluons 0,[6|. Technically this can be formulated by introducing 
soft as well as ultrasoft gluons in the theory 0]. 

In Refs. js], H, [lo|, based on the hierarchy m ^ mv ^ mv'^, it was suggested to employ a 
series of EFT's, one for the scales m > fi > mv and one for mv > fi > mv"^, called pNRQCD 
("potential" NRQCD). The pNRQCD is derived in a two-step matching procedure where the 
intermediate soft matching scale and the pNRQCD renormalization scale are independent 
at the field theoretic level. This is appropriate e.g. for static quarks where the momentum 
is fixed by the quark separation r ~ 1/mv and not correlated with the dynamical energy 
fluctuations E ~ mv'^. The potential interactions are generated from integrating out soft 
fluctuations at the soft matching scale. 



In Ref. is was pointed out that for dynamical heavy quarkonium systems the dis- 
persive correlation between ultrasoft energy and soft momentum scales, E ~ p'^/m, can be 
implemented systematically at the field theoretic level by matching to the proper EFT di- 
rectly at the hard scale /i ~ m (one-step matching). The EFT, called vNRQCD ("velocity" 
NRQCD) has a strict power counting in v. It contains soft and ultrasoft degrees of freedom 
as well as soft and ultrasoft renormalization scales, and fiu- Through the different ve- 
locity counting of soft and ultrasoft fields they are correlated as fiu oc where usually 
the constant of proportionality is set to 1, i.e. fiu = A'-l/''^ = mu'^. The renormalization 
group running of the theory is then conveniently expressed in terms of the dimensionless 
parameter u. The matching is carried out at the hard scale v = 1 and the theory is evolved 
to ~ f ~ cts of order of the relative velocity of the two-body system for computations of 
matrix elements. The runnin g pr operly sums logarithms of the momentum and the energy 
scale at the same time [nl, [HJ Il3l . Il4| and is referred to as the velocity renormalization group 



(VRG) [llj. Within dimensional regularization the powers of /x^ and /i^ multiplying the 
operators of the Lagrangian are uniquely determined by the v counting and the dimension 
of the operators in d = 4 — 2e dimensions. 

In Ref. p^], where a three-loop renormalization of quark-antiquark vertex diagrams was 
carried out for fermion-antifermion S-wave states, is was demonstrated that the one-step 
matching principle, upon which vNRQCD is based, is consistent under renormalization at the 
subleading order level, when sub divergences need to be subtracted and when UV divergences 
from heavy quark-antiquark loops and from soft or ultrasoft gluons appear simultaneously in 
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single diagrams. Up to now there is no analogous demonstration for the two-step matching 
principle. 

With the VRG, the running of potentials relevant for the next-to-next-to-leading logarith- 
mic (NNLL) description of the nonrelativistic dynamics of pairs of heavy quarks and colored 
scalars has been determined in Refs. H H Q [13 and [l9|. The NIL running of lead- 
ing (in the v expansion) S-wave currents for heavy quark- ant iquark production and of the 



leading S- and P-wave currents for pairs of colored scalars were determined in Refs. |14l . Il8 



and |19|]. Corresponding work in the pNRQCD formalism for heavy quarks was carried out 
in Refs. (20I . [21I . 22, 23|. In these computations, diagrams which simultaneously contain 
both UV divergences from heavy quark-antiquark loops and from soft or ultrasoft gluon 
loops do not arise. The final results based on the VRG and on pNRQCD agree, after the 
scale correlations from the one-step matching procedure of vNRQCD are imposed in pN- 
RQCD. Concerning the NNLL running of the production currents, only the contributions 
from three-loop vertex diagrams for the _leading currents describing quark-antiquark pair 
production in an S-wave configuration 
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are known at this time in vNRQCD. For the 
NNLL contributions that arise from the subleading order evolution of the coefficients that 
appears in the NLL anomalous dimension of the current at present only the results for the 
Coulomb potential [13, lO], the spin-dependent 1/m? potentials [l^, 25 1 and from ultrasoft 
quark loops [26| have been determined. 

The main phenomenological application of the nonrelativistic EFT for the situation 
m ^ mv 3> mf ^ > Aqcd is top quark pair production close to threshold at a future 
Linear Collider [23]. Here the summation of logarithms of v is crucial to control the large 



normalization uncertainties |28l . |29j that arise in fixed-order nonrelativistic computations 



that only account for the summation of Coulomb potential insertions |30|, l3l|]- One can ex- 
pect that the summation of logarithms of the velocity is also important for the description 
of squark pair production at threshold 19|]. Recently, it was also found that it is crucial to 
apply the EFT for predictions of the e^e~ ttH cross section at ^/s = 500 GeV [H, 3^ 
since here the tt final state interactions are dominated by nonrelativistic dynamics. 

In this work we address the form and construction of non-relativistic interpolating cur- 
rents for production and annihilation of a color singlet heavy particle-antiparticle pair with 
general quantum numbers '^^'^^Lj in n = d— 1 = 3 — 2e dimensions for quarks and colored 
scalars within NRQCD. In general, the interpolating currents for a specific "^^^^Lj state are 
associated to irreducible tensor representations of the SO(n) rotation group that can be built 
from the particle-antiparticle relative momentum and the bilinear covariants of the particle- 
antiparticle field operators. While in three dimensions the basis of interpolating currents 
for the different '^^^^Lj states is, by construction, unique, in general n 7^ 3 dimensions for 
fermions there exist evanescent operator structures that make the choice of basis for the 
interpolating currents ambiguous. This is because the structure of irreducible SO(n) repre- 
sentations for n 7^ 3 is inherently different and in general richer than in three dimensions. 
Technically the ambiguity is related to the number of sigma matrices that generate the SO(n) 
rotations for the spinor wave functions. For explicit computations a specific choice of basis 
has to be made, which is associated to a specific choice of a renormalization scheme. A very 
similar problem has been treated systematically already some time ago in the framework of 
relativistic QCD and the effective weak Hamiltonian [34, 35, 3^, and much of the statements 
that apply for the relativistic case can be directly transferred to the nonrelativistic theory. 
However, due to the nonrelativistic power counting some even more specific statements, 
can be made. In particular, the NLL order matching conditions and anomalous dimensions 
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of the ^'^"'"^Lj interpolating currents that are leading in the nonrelativistic expansion are 
scheme independent. 

Using the explicit form of the interpolating currents obtained in this work we also deter- 
mine the NLL anomalous dimensions of the leading (in the v expansion) currents describing 
color singlet heavy quark-antiquark and squark-antisquark pair production for arbitrary 
spin and angular momentum "^^^^Lj configurations. The results for low angular momentum 
states (S,P,D) are e.g. relevant for angular distributions at the threshold or for production 
and decay rates of squark-antisquark resonances in certain super symmetric scenarios where 
squarks can have a very long lifetime. The results also shed some light on the importance of 
summing logarithms of v for the production and annihilation rates of high angular momen- 
tum states. In this respect our results help to complete analogous higher order considerations 
in previous literature for the energy levels [37|] and the wave-functions [2^ 



For the presentations in this work we employ the notations from vNRQCD based on a 
label formalism for soft fields and the heavy quarks (or scalars) We note, however, 

that most of the results obtained in this work are applicable to NRQCD in general. In 
particular, the NLL order anomalous dimensions obtained here are also valid in pNRQCD 
once the vNRQCD scale correlations from the one-step matching are imposed. 

The outline of the paper is as follows: In Sec. [Ill we discuss the form of the spherical 
harmonics in n dimensions and present the form of currents describing states with arbitrary 
angular momentum L and total spin zero. We also present the n-dimensional form of the 
Legendre polynomials and demonstrate in an example why the n-dimensional form of the 
spherical harmonics is needed in dimensional regularization. In Sec. IIIII we discuss the 
properties of cr-matrices for n 7^ 3 and the importance of evanescent operator structures 
that can be built from the cr-matrices. We present a simple basis of interpolating currents 
describing fermion-antifermion pairs in a spin triplet state and having arbitrary L and total 
angular momentum J. In Sec. |IV]we determine the anomalous dimension of the current for 
spin singlets, and in Sec. |V] those of the current for spin triplets. Section |V] also contains 
a discussion on the scheme-dependence of the spin-dependent potentials. In Sec. |Vl] we 
determine and solve the resulting anomalous dimensions for heavy quarks and colored scalars. 
The results are analysed numerically for a few cases. In Sec. I VIII we comment on the scheme- 
dependence of results that can be found in previous literature. The conclusions are given 
in Sec. IVIIII We have added a few appendices: Appendices \M and [B] contain the derivation 
of the tree level NRQCD matching conditions for the well known processes qq —>■ 27 and 
37 accounting properly for the existence of evenescent operator structures. In App. [Cisome 
more details are given on the determination of the spin triplet currents shown in Sec. |Vl We 
also derive an alternative set of currents that is equivalent for n = 3 but inequivalent for 
n ^ 3. Finally in App. [D]we give results for the UV-divergences for the elementary integrals 
that are needed for the determination of the anomalous dimensions of the currents. 



II. SPIN SINGLET CURRENTS 



In this section we discuss the form of the interpolating currents describing the production 
of a particle-antiparticle pair with total spin zero for arbitrary relative angular momentum 
L (^^~^^Lj = ^Ll). They are relevant for pairs of scalars or for quark-antiquark pairs in a 
spin singlet state. The currents naturally have their simplest form in the cm. frame where 
relevant angular dependence can only arise from the cm. spatial momenta of the pair. Thus 
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the generic structure of the production currents is 

ijl{x)r{p)x%{x), (1) 

where r(p) represents an arbitrary tensor depending on the cm. momentum label p. For the 
case of scalars ipp and Xp ^i-re scalar fields with ip^ = V^p, while for fermions ipp and Xp ^i-re 
Pauli spinor fields with x*-p = (^cr2)X-p- (Note that we adopt the standard nonrelativistic 
spinor convention of antiparticles with positive energies. In this convention fermions and 
antifermions have the same spin operators.) The corresponding annihilation current can 
be obtained by hermitian conjugation. The interpolating currents associated to a definite 
angular momentum state L are related to irreducible representations of the tensor F with 
respect to the rotation group SO(n). Since the transformation properties of the fields are 
not relevant in this respect for the spin zero case, it is sufficient to identify the irreducible 
tensors that can be built from the spatial momentum vector p^, i = 1, . . . ,n, where n = d — 1. 
The irreducible tensors associated to the angular momentum L are up to normalization just 
the spherical harmonics of degree L. In the following we give a general discussion in n 
dimensions. All results reduce to the well known results for n = 3. 



Spherical Harmonics 

The spherical harmonics of degree L are the polynomials u{x.) in M" which are homoge- 
neous of degree L (i.e. u{rx.) = r^u{x.)), harmonic (i.e. they satisfy the Laplace equation 
A]Knu(x) = V^m(x) = 0) and restricted to the unit sphere 5'"'"^. The spherical harmonics 
Y;,A/(n, x) of degree L, with M = 1, . . . ,ni and x = x/|x|, form an orthogonal basis of a 
n^-dimensional vector space with 

fn + L-l\ fn + L-3\ T{n + L-2) 

"^=1 L )-[ L-2 ) =('^ + "-^) F(n-l)F(L + l) - 

Any different iable function on the unit sphere S"""^ in can be expanded in terms of a 
convergent series of spherical harmonics. A representation in terms of cartesian coordinates ^ 
of the spherical harmonics of degree L is given by the totally symmetric and traceless tensors 



with L indices T*i'"*^(x), where the indices ii,...,iL are cartesian coordinates [40]. The 
restriction to the unit sphere is not essential regarding the transformation properties under 
SO(n) rotations and can be dropped for our purposes. An explicit expression for T*^ ' *^(a;) 



^ Although it is possible to write down results with spherical coordinates for integer n (see e.g. Ref. [39|), 
it is more convenient in practice and for actual computations to use cartesian coordinates. 
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reads: 



X . . .X 



(2L + n-4) (2L + n-4)(2L + n-6) 

. + (_1)[V2]. 



(x2)i^/2]$;--(x) 



(2L + n - 4) . . . (2L + - 2 - 2[L/2]) 



[L/2] 
fc=0 

. ^^ i.r(^ + L-A;-r 



r(f + L-i) ' 



incqu. pcrmut. ^ 
(pi...Pi)of 



= 5'^'^ . . . 5^2fc-ii2fc ^i2k+i __^iL ^ (inequ. permut.'s) , (3) 

where the symbol [. . .] denotes the Gauss bracket. The sum in the tensors extends over 
all sets of k pairs of indices that can be constructed from L indices (L > 2k)] the number of 
terms in the sum is thus (^) (2k — 1)!!. Note that we call two permutations equivalent if they 
lead to the same term in the sum. Although we believe that the expression for T*i'"*^(a;) 
has been shown somewhere in the literature before, we were not able to locate a suitable 
reference. A few simple examples are 

x^ 

r(x) = x\ r^(x) = x'x' — 

n 

x2 

T^^ipi) = x'x^x^ {x'5^^ + x^5''' + x''5'^) . (4) 

n + 2 

From the Laplacian 



c?x*^ \dr'^ ^ r drj 

_L2 = —x''—x^ - —x''—x^ = x^V^ - xV^ — -(n- , (5) 

dx^ dx^ dx^ dx^ dx^ dx^ dx^ ' 

one can find an explicit form for the squared angular momentum operator and, using the 
homogeneity relation AjRnT*^'' *^ = 0, the eigenvalue equation 

L2 = L{L + n-2) T'-'^{x.) . (6) 

We define the interpolating currents that describe production of a spin singlet and angular 
momentum L state {^^~^^Lj = ^Li) as 

(jf=°)"-^" = ^l{x)T^-^-{jp)x%{x). (7) 

The current in Eq. ([7]) is the dominant ^Ll current in the f -expansion. Higher order currents 
are obtained by including additional powers of the scalar term p^. 
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As an example, let us consider the currents with angular momenta S, P and D. They 
are relevant in the electromagnetic production of heavy charged scalars from e^e~ and 77 
collisions. The e+e~ annihilation at lowest order in the expansion of the electromagnetic 
coupling produces P-wave states: 



iMe 



^-i-^Q,v{k')^'u{k)T\^), 



(8) 



with /c, k' the momenta of the incoming e+e~ in the cm. frame and Qq denoting the electro- 
magnetic charge of the scalars. Note that potential color indices are always implied. Pairs 
of heavy scalars in S- and D-waves are first produced in 77 collisions. The first terms of 
the amplitude in the expansion in the cm. three-momentum of the outgoing particles read: 



iM-yy = 32i TTCK 



4 2w? 



= -8i7rQ;g^eie^2 1 



2p^ 



m 



2 Q^i47^^^(p) + ---- 



(9) 



The first term in the last equality is the S'-wave contribution (T = 1) whereas the second is 
the D-wave contribution described by the spherical harmonic of degree 2, T'^'2(p). 
Some useful relations for the T tensors read: 



L + n-3 



2L + n - 4 



.2 L 

(2L + n-2)(2L + n-4) ^ 



x) 



3< 
J = l 



2L-Fn - 2 ^ 



x), 



d_ 



^i^-"(x) = $^<5^^^$r"^-'^(x) = $--^(x)-x^$-^-(x) 



i=i 

where we use the notation that hatted indices are dropped, and A^'^'^ = ^{A'^ - A^'). 



(10) 
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Legendre Polynomials 

The spherical harmonics of degree L satisfy the addition theorem 



riL 



V YLMin, x) Yl^n, y) = ^ Pi(n, x.y) , (11) 



a. 

M=l 



where PL{n,t), t G [—1,1], is the Legendre polynomial of degree L generalized to n di- 

r(n/2) 



mensions, with PL(ra, 1) = 1 and cr„ = yt~M ^^^^ ^'^^ ^^^^ sphere S*" ^. The 



Legendre polynomial is related to the Gegenbauer polynomial C^(t).^ One has PL{n,t) = 
^^rn^J^n-'^ ^L (^) i where A = The exphcit form of the generalized Legendre polynomial 
reads Ell, 



P(nf)- r(n-2)r(L + l) {-if T [r^ + L - k) 

^L[n, t) - ^^^^ r(n + L - 2) ^ r(A; + 1) r(L - 2A; + 1) ' ^'^^ 

From Eq. (fTTj) follows an addition theorem for the T tensors upon total contraction of their 
indices: 

Ti,„(x,y) = T^i-*^(x)T^i-*^(y) = |x|^ |y|^ iV^ Pi(n, t) , = ^\ ' (^3) 

where 

iV,^ r(L + n-2)r(^) 
^ 2^r(L + f -l)r(n-2) ^ ^ 

Another very useful contraction is: 



Consistency Consideration 

The use of the generalized currents in Eqs. Q and ([7]) is mandatory to obtain con- 
sistent results in dimensional regularization in accordance with SO(n) rotational invari- 
ance. As an example let us consider the nonrelativistic three-loop vacuum polarization 
diagram shown in Fig. [1] with two insertions of the Coulomb potential. Inserting the cur- 
rents that produce and annihilate the particle-antiparticle pair with angular momentum L 
and fully contracting the indices, the quantity we want to compute is, up to a global factor 

Nl J (qf + 5) (qi - q2)2 (q^ + 5) (qa - qs)^ (q^ + 5) 



^ The index A is chosen such that the are orthogonal on the n-dimensional sphere 
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FIG. 1: Three-loop diagram with two insertions of the Coulomb potential (black dots). 



with 6 = —rriE — ie. We can now proceed in two different ways to do the computation. From 
Eq. (JT3ll . the contraction of the T's at the ends generates a Legendre polynomial depending 
on the angle between the loop momenta on the sides: 

A['^ = [ D-a, D"q, D-as . . .... • (17) 



(q? + 5) (qi - q2)2 (qi + 5) (q. - qs)^ (qi + 5) ' 

On the other hand, we can first shift the loop momenta dependence of the current T*^" *^(qi) 
by using the relation 

j (i'^qiT^- -^^(qi)/(qi,q2) = AT^--^-(q2), 

A = y"ci^qi|qi|^|q2r''PL(ri,qi.q2)/(qi,q2), (18) 

which is a consequence of rotational invariance if /(qi, q2) is a scalar function depending on 
q?) qi) qi ■ qs- This gives an alternative expression with two Legendre polynomials: 

/!"■>- ^D-n D"n D'n llA'-M'' Pdn, <li.<h) P^n, qiAs) , „, 



Had we worked in n = 3 from the beginning, we would have written down Eqs. (ITTI) and f[T9|) 
with the generalized Legendre polynomials replaced by their n = 3 expressions. Let us call 
the corresponding expressions as 2 . Since the integrals are power divergent, we 

can compute them using dimensional regularization in = 3 — 2e dimensions and check 
whether they give the same result. The first non-trivial case is L = 2, because one has 
Po,i(3,a;) = PQ^i{n,x). The results for ^4^ , A^2 then read 

4'^ = J^W^-Vr{l5 + 47r^ + 6[327r^-120C(3) + 23l]+O(6^)}, 

Af = J^(5)^-Vr{l8 + 47r2 + 6[327r2-120C(3) + 264-41og8]+O(e2)}. 

(20) 

A difference arises even in the first term in the e expansion, which reflects the fact that 
2 are not the same quantity. The reason is that the shift (fT8|) performed in n = 3 
dimensions and the following evaluation using dimensional regularization do not commute 
if the integrals are divergent. Rotational invariance in n dimensions must be satisfied at 
every step if dimensional regularization is used. This is automatically achieved by using the 
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SO(n) tensors T as shown in Eqs. (fT7|) and (fT9|) . The computation with generahzed Legendre 
polynomials corresponding to A^^'^ and A^^'^ in Eqs. fll7lll9p can be shown to produce the 
correct outcome. For L = 2 the correct result reads 

■^^'^ = ^^'^ = 3^ ('^)^~'/^^^ {l5 + 47r2 + e [32n' - 120 C(3) + 236] + 0{e')} , 

which disagrees with A\ in terms proportional to e"^ for m > 0, and with A2 for all terms 
in the e expansion. 



III. SPIN TRIPLET CURRENTS 

In this section we discuss the form of the interpolating currents describing the production 
of a fermion-antifermion pair in a spin triplet S = 1 state for arbitrary relative angular 
momentum L (^^~^^Lj = ^Lj). As in Sec. [TTlthe currents are defined in the cm. frame. 



Pauli Matrices 



Since the treatment of Pauli cx-matrices in n dimensions involves a number of subtleties, 
we briefly review some of their properties relevant for the formulation of the currents. Many 
of the properties can be directly obtained from the corresponding properties of the 7-matrices 
in d = n + 1 dimensions. The a-matrices a* {i = 1,. . . ,n) are the generators of SO(n) 
rotations for spin 1/2. They are traceless, hermitian, satisfy the Euclidean Clifford algebra 
{(J*, a^} = 26^\ and can be defined such that cr^ = cr* = —a^cro^- While traces of a product 
of an even number of cr-matrices in n dimensions can be expressed as a multiple of Tr[l] using 
the anticommutator, traces of a product of an odd number are identically zero and for the 
case of three cr-matrices can require additional rules to yield results that are consistent with 
the relations known from n = 3. In this respect the product of three different cr-matrices 
can be somewhat considered the three-dimensional analog of 75 in four dimensions.'^ As for 
the case of the 7-matrices ^] products of cr-matrices in arbitrary number of dimensions 
cannot be reduced to a finite basis, but represent an infinite set of independent structures. 



In analogy to Ref. [3J] it is convenient to define 

^n-i™ ^ (m = l,2,...) (21) 



^ The vanishing of traces of a product of an odd number of cr-niatrices in n dimensions is related to 
the simultaneous use of the cyclicity property of the trace operators and the Euclidean Clifford algebra 
relations for all cr-matrices. Products of an odd number of cr-matrices, however, do not arise from insertions 
of potentials which we consider in this work. Such traces can, however, arise e.g. when the spin-dependent 
radiation of ultrasoft gluons or annihilation potentials need to be considered. Thus, in general, the cyclicity 
property of the trace operation should only be used after renormalization is completed - if the Euclidean 
Clifford algebra is applied for all cr-matrices. 
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for the averaged antisymmetrized product of cr-matrices, so e.g. a*-' = l/2((jV-' — a-^o"*). It 
is straightforward to derive the following useful relations: 



m 

-I)'" CT^^i-^™ + ^(-1)^' 6^'^ a^^ -^^-*™) , (22) 



= 4(-l)'™ ^ (-1)'=+^ (^^'^^a^^---^*---^™^ - (5^*'=a^i---^'=---^'"^') , (24) 



A;=l 



- (aV^i-^'-a^ + aV^-^-a^) = (-I)'" 5^'^ 
2 



A:=l 

^ ffV^i -^-a^' = (-1)" (n - 2m) cx*^-^"' , (26) 

A: 

a'^-"^- a*^-"^ = (-1) 2 ri...(n-A; + l) = (-1)^ ^ ^ ' 



T(n-k + l)' 
ii...tk 

(27) 

Tr (cT^i ...CT*2™) = Tr[l] ^ 5'piip2 . . . ^W-iW . (28) 

incqu, pcrmut. 
(Pl"-P2m) °f (l--.2m) 

In the last relation, 5p is the signature (sign) of the respective permutation {p_^p^" ), and 
for each 5*™*" only m < is allowed. Note that we call two permutations equivalent if they 
lead to the same term in the sum. For n = 3 we have a*-' = ie^^^a^, a^^^ = ie^^^l, and 
^«i - «m _ g fQj, m > 3. Thus the latter are evanescent for n ^ 3. For m < 3 the cr*^ "*™ are 
related to physical spin operators. 



S-Wave Currents 

The general interpolating current describing the production of a fermion-antifermion pair 
in an S-wave state and in an arbitrary spin state in n dimensions has the form 

i3L=oY''-'-^ = ^lix) a'--'- i^cr2)X-pix) . (29) 

A SO(n) rotation by an angle 6 around the axis n leads to the transformed current 
Dei'^la''-'^" (z(T2)xlp) = ^J,5(^)t(T^i-*™ S{9){ia2)x-p, where S{9) = exp(-z^n.o-/2). Since 
5 (6') 5(0) = D^^{6)(7^ , where D^^{6) is the rotation matrix for n- vectors, the currents in 
Eq. fl29|) are tensors with respect to SO(n) rotations. The tensors are irreducible due to the 



antisymmetry of the cr*^ "*™ [40| and each have (J^) independent free components. Their eigen- 
values with respect to the square of the total spin operator = [a-p lap + Ip (8> (Tapf/A — 
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{n[la ® lap] + [cTp ® crap])/2, wheie the indices of the cr-matrices are summed over, read 



= ^ (n + (-1)" (2m - n)) (jL=o)f^^ -*'"^ • (30) 

For the physical currents with m = (0,1,2,3) one thus finds the spin eigenvalues {0,n — 
1, 2, n— 3). Note that the spin eigenvalue for the unphysical S-wave currents a*!" '™ for m > 3 
are non-zero. While for n = 3 the spin singlet currents for m = 0, 3 are equivalent, and 
likewise the triplet currents for m = 1,2, each of the currents represents a different irreducible 
representation of SO(n) for n ^ 3. It is a very instructive fact that the action of onto 
the singlet current j]*L*o*^^ does not give zero for n ^ 3. Thus to achieve that spin-dependent 
potentials do not contribute for physical predictions involving spin singlet currents, in general 
additional finite renormalizations are required in analogy to Ref. ^3J], unless a scheme for 
potentials or currents is chosen such that they vanish automatically. However, even if such 
a scheme is adopted, matrix elements, matching conditions and anomalous dimensions can 
depend on the spin-dependent potentials at nontrivial subleading order^. 

All of the physical currents can arise in important processes. In Tab. [T] the leading or- 
der matching for a number of different currents is displayed. In general there are several 
relativistic currents that can lead to the same nonrelativistic current [42|. Note that the 
respective production and annihilation currents are related either by hermitian or antiher- 
mitian conjugation. Except for 75, which is treated as fully anticommuting (75 = (JJ)), three 
dimensional relations have not been used. The full theory current with the 7-structure Y''^ 
for example arises for fermion pair production in 77 collisions, as shown in Appendix 1X1 
However, also evanescent currents naturally occur in standard processes, such as the current 
Jl=o*^' that arises in fermion- ant if ermion pair annihilation into three photons. The explicit 
computation can be found in Appendix [Bl Also note that the differences between the two 
different singlet and triplet currents correspond to evanescent operators as well. It is well 
known from subleading order computations based on the effective weak Hamiltonian that 
one needs to consistently account for the evanescent operator structures that arise in matrix 



elements of physical operators when being dressed with gluons. In Ref. [3J] it was shown 



that a renormalization scheme can be adopted such that a mixing of evanescent operators 



into physical ones does not arise. Moreover it is also known [35|, |36| that modifications of 
the evanescent operator basis (e.g. by adding physical operators multiplied by functions of e 
that vanish for e — > 0) and similar modifications to the physical operator basis correspond 
to a change of the renormalization scheme. While this does not affect physical predictions, 
it does affect matrix elements, matching conditions and anomalous dimensions at nontrivial 
subleading order. Thus precise definitions of the schemes being employed have to be given 
to render such intermediate results useful. 

In the framework of the nonrelativistic EFT these properties still apply. However, the 
velocity power counting in the EFT allows for even more specific statements. Concerning 
interactions through potentials, transitions between S-wave currents in Eq. fl29|) that are in- 
equivalent cannot occur because the potentials are SO(n) scalars and the currents are (due to 



^ To be more specific, we refer to orders of perturbation theory where subdivergences need to be renormal- 
ized. 
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Full theory T 


r(p,cr) 


/ 


^ ' m 


ji...i2fc+i 




7 7 


W I ^) 2^.7=1 1 m 


^0^n...i2fc+i 


(★) (_i)fecrn-*2fc+i 




(*) (_i)fccjn-i2fe 


^jl---«2fe+1^5 




^0^n..i2fc^5 




^0yi...i2fc+i^5 





TABLE I: Nonrelativistic production currents ^pr(p,<T) [ia2)X-p arising from the leading order 
matching with the full theory currents ipTiJj. The nonrelativistic normalization of the full theory 
spinors has been used (see Eq. (|A1|1 ). The notation stands for the averaged antisymmetrized 

product of 7-matrices for indices other than zero. For the corresponding annihilation currents 
X^p(— icr2) r'(p, cr) -(/jp one has F' = — F for all entries with the (*) symbol and F' = F for the 
others. 



the different symmetry patterns of their indices 40|) inequivalent irreducible representations 



of SO(n). Even for currents with L ^ and for the spin-dependent spin-orbit and tensor 
potentials (which is all we need to consider at NNLL order) one can show with Eqs. flMj) 
and fl2Sl) that transitions between currents containing cr*^"'*™ with a different number of 
indices cannot occur (see Sec.|Vll). Transitions can, however, arise between currents with 
different angular momentum, such as for the tensor potential that can generate transitions 
between L and L' = L ±2 (see Sec.|V]). The same arguments apply to the exchange of soft 
gluons (in the framework of vNRQCD) since they cannot appear as external particles and 
furthermore need to be exchanged in pairs due to energy conservation. In this respect the 
effects from soft gluon exchange effectively represent modifications to the potential interac- 
tions (see e.g. Ref. [lH). Concerning the exchange of ultrasoft gluons, transitions between 
currents containing cr*! ' *™ with a different number of indices can arise, but only if the in- 
teraction is spin-dependent. The dominant among these interactions corresponds to the 
operator ip^a^^ y ippA^ , where is an ultrasoft gluon. This operator can only contribute at 
N^LL order, which is beyond the present need and technical capabilities. Thus in practice, 
transitions between currents containing cr*^ "*™ with a different number of indices do not 
need to be considered. 

So for the S-wave currents in n dimensions one can employ either one of the two spin 
singlet [k = 0,3) or triplet currents {k = 1,2) in the EFT and the difference corresponds 
to a change in the renormalization scheme. This means in particular that as long as the 
renormalization process is restricted e.g. to time-ordered products of the currents, one can 
(but does not have to) freely use three-dimensional relations to reduce the basis of the 
physical currents. However, once the basis of the physical currents is fixed (where each 
current is irreducible with respect to SO(n)), one has to consistently apply the computational 
rules in n dimensions discussed in the previous sections. As an example, instead of using the 
current ip^a'^^^{ia2)X-p that arises in 77 collision, one can employ the current ipp{'i<^2)x*-py 
defined in n dimensions, times the e-tensor ie^^'' defined in three dimensions, which means 
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that the e-tensor is zero if any of its indices takes a value different from 1, 2, 3. Concerning the 
e-tensor, this works because the e-tensor does not play any role during the renormalization 
procedure as long as we only consider time-ordered products of the currents.^ This justifies 
the approach in Sec. [IT] where we have only considered spin singlet currents for fermions 
involving the current ?/'j,(icT2)xlp- This scheme is also advantageous practically because 
for the current 'ippO'^-''' {icr2)x*-p the matrix elements of the spin-dependent potentials that 
vanish for n = 3 can, as discussed above, give evanescent contributions for n ^ 3. Moreover, 
from these considerations we can also conclude that currents containing evanescent cr*^'"*"' 
matrices (m > 3) can (but do not have to) be dropped from the very beginning in the EFT 
as long as one does not need to account for spin-dependent ultrasoft gluon interactions. 

An important lesson to learn from this discussion is that partial results at nontrivial 
subleading order obtained from the threshold expansion j43| of full theory diagrams such 



are 



as e.g. the contributions from the hard regions obtained in Refs. |44J . I45l . 1461 . |47[, 
equal to the EFT matching conditions only in schemes for the effective theory currents and 
potentials that are compatible with the nonrelativistic reduction of the 7 matrices that has 
been used during the threshold expansion. In general, there are additional contributions to 
EFT matching conditions to account for the scheme choices made in the EFT. Note that in 
the threshold expansion different scheme choices can also be possible, in particular for the 
treatment of 75. In Sec. IVIII we comment on scheme dependences of a number of partial 
results that can be found in the literature. 

At this point one might also ask whether 75 needs special treatment in the nonrelativistic 
EFT in n dimensions. Chirality and the flavor symmetries are not relevant in the EFT for 
a single heavy particle-antiparticle system and their potential effects are contained in the 
matching contributions of the EFT to the full theory. For the matching relations shown in 
Tab. [Hwe have used a totally anticommuting 75. Since the resulting effective theory currents 
have well defined SO(n) transformation properties and reduce to the proper nonrelativistic 
currents for n ^ 3, this represents a consistent scheme choice from the point of view of the 
effective theory computations. A different ansatz for 75 such as 75 = i'-f^'j^'-f'^'j^ (which is the 
fully consistent one in the full theory) just corresponds to a different choice of scheme for 
the currents in the nonrelativistic EFT, and both schemes can be used consistently. Note 
that from the results in Tab. [I] without an explicit 75 one can also derive the form of the 
EFT currents for the ansatz 75 = ij^j^j'^'j^ (see also Sec. IVIip . 



Arbitrary Angular Momentum 

Based on the results obtained in the previous sections it is straightforward to construct 
spin-triplet currents with arbitrary angular momentum L {^^^^Lj = ^Lj)). They can be 
obtained by determining irreducible SO(n) representations from products of the tensors 
T*i "*^(p) describing angular momentum L and the spin-triplet S = 1 currents discussed in 
the previous paragraphs. Due to the different symmetry patterns of the symmetric T*i'- *^(p) 
tensors and the antisymmetric 'j^cr*i" **(i(T2)xlp currents one needs to apply the construction 
principles for general tensors which state that tensors are irreducible with respect to 



^ This procedure cannot be applied in this simple form, if the initial state that is involved in the quark- 
antiquark production process is also involved in the renormalization procedure, as can be the case for 
QED corrections. 
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SO(n) if and only if they are traceless and if their indices have a symmetry pattern according 
to a standard Young tableau. As for the case of the S-wave currents the physical basis for 
arbitrary spatial angular momentum is not unique due to the existence of evanescent operator 
structures. 

Here we construct currents with fully symmetric indices because of their comparatively 
simple form and because their number of indices is equal to the total angular momentum 
J quantum number. For the case J = L ± 1, these currents are contained in the reduction 
of the reducible tensor A'^'^""^^+^ = '?/'^T*^" *^(p)cr*^+i(i(T2)x-p- Upon symmetrization and 
removal of all traces one obtains the '^^'^^Lj = ^L^^i current 



L+1 L+l 

rf = i(p,<T)*l-^^ + l = T*l-^'=-*^ + ^(p)(T^'= §ikijrpil-t,...T,...il^ + li(^p^^i _ 

^ J-j I 2 
fc=l k<J 

fe=l 

(31) 

The subtractions in the second term of Eq. fIST]) needed to achieve zero traces are for them- 
selves "^^^^Lj = ^Ll-i currents, so we can define 



ri=Up,cr)*i-"^-i = T''-'^-'\p)a\ (32) 

The tensor ' also contains ^Ll currents. However, their indices are not fully symmet- 
ric. Their form is discussed in the Appendix IB] together with a more detailed derivation of 
the two currents displayed above. A fully symmetric ^L^ current can be obtained from the 
reduction of the tensor i?*i---*i'+2 = -j/jj^ T*i - *^ (p)(T*-^+i*^+2(i(T2)xlp. Upon removal of traces of 
j^H---tL+2 Qj2g QQj^ identify the ^L^ current 



rf=^(p,o-)^i-^^ = ^ T*i-^'=-^^'^(p)a^^'= . (33) 
fe=i 

The tensor _B*i - *i+2 also contains evanescent currents and also currents describing ^Ll±i 
states that differ from the currents given above for n 7^ 3. Their form is also discussed in the 
Appendix [Bl The number of independent components of the spin-triplet currents defined 
above is that of a fully symmetric tensor with J indices, i.e. equal to nj given in Eq. ([2]). 
All three currents have parity and charge conjugation P = (— l)"^"*"^ and C = (— 1)^"'"'^, 
respectively. Also note that that r£^J(p, cr) and ir|^^(p,cr) are hermitian matrices. 

The eigenvalue equation of the currents with respect to the spin-orbit operator S.L = 
-[S^, S''] pid/dp^ = -l/4( [al, a^] ® lap + Ip ® [cr^^, o-^J ) pid/dp^ has the form 



= -\ {[^^ r(p, <Tr-^{ia,) + A r(p, ar-^ {^a,)[<T\ ^T)x-p , (34) 
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„-5=0 
Jj=L 


.■s=i 


^S=l 
Jj=L 


.■5=1 
JJ=L+1 




L{L + n-2) 


L(L + n-2) 


L(L + n-2) 


L{L + n-2) 







n — 1 


2 


n — 1 


S.L 





-(L + n- 2) 


-(n-2) 


L 


J2 


L(L + n-2) 


J{J + n-2) 


J(J + n - 2) + 6 - 2n 


J{J + n- 2) 



TABLE II: Spin and orbital angular momentum quantum numbers of the currents defined in 
Eqs. (j7l3f I32I33P in n dimensions. 



where r(p, cr)*! "*-^ represents any of the tensors between the fermion fields. The eigenvalues 
of the different currents defined in Eqs. (17|31|32II33I) with respect to the operators L^, S^, 
L.S and are summarized in Tab. [TTl For the computation of the NLL order anomalous 
dimensions of the currents from time-ordered products of two currents (Sees. [IV] and |V]), the 
total contractions of spatial and spin indices (i.e. taking the trace of products of a-matrices) 
are useful: 



Tr 



r^=}^(p, (7)''-^^+' r^=y^ „\n-iL+, 



Tr 



r 



5=1 1, 



L+iiq,'^j 

^5=1, 



Tr 



5=1/ 



(L + 1) (2L + n) TlAv. q) Tr [1] , 

L{L + n- 2) Ti,„(p, q) Tr [1] , 
Ti,„(p,q)Tr[l], (35) 



where the function Tl^„ is defined in Eq. ( !T3|) . Note that the expressions only involve traces 
of products of an even number of cr-matrices which can be computed unambiguously from 
the anticommutator of two a-matrices. 



IV. NLL ORDER ANOMALOUS DIMENSIONS: SINGLET CASE 

For the computation of the renormalization constants of the physical currents defined in 
the previous sections one has to determine the overall UV- divergence of particle-antiparticle- 
to-vacuum on-shell matrix elements of the currents from vertex diagrams. However, this 
method is complicated by the fact that for on-shell external particles with E = p^/m the ver- 
tex diagrams also contain IR-divergent Coulomb phases that have to be distinguished from 
the UV-divergences. To avoid this issue we use the method from Ref. [l5| which considers 
current-current correlator graphs obtained from closing the external particles-antiparticle 
line of the vertex diagrams with an additional insertion of the (hermitian conjugated) cur- 
rent. This means that one has to determine diagrams with one more loop, but in this way 
all IR-divergences associated to the Coulomb phase cancel and one needs to compute a fewer 
number of diagrams. Technically, the UV-divergences associated to the currents appear as 
sub divergences of the correlator diagrams, which can, however, easily identified from the 
absorptive parts of the correlators. 

Since there are no UV divergences in the EFT that can lead to a LL anomalous dimension 



of the currents that are leading in the nonrelativistic expansion [ll|], the NLL order level 
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ocx> 

(d) (e) 5ci 

(^^^^^^^^^^ 

FIG. 2: Three-loop current correlator diagrams and the counterterm diagram for the computation 
of the NLL anomalous dimension of the currents. Combinatorial factors are suppressed. 



represents the lowest order in which these currents are renormalized. The corresponding 
three-loop current correlator diagrams are depicted in Fig.[2l The terms vjf"^ and Vg*'', Vr^^ 
are the color singlet Wilson coefficients of the spin-independent potentials of order asV~^ 



(1/k^) and a^f (1/m^, (p^ + p'^)/(2m^k^)), respectively [l6|,ll7(], while are the Color 

singlet Wilson coefficients of the sum operators C^*'', O^j}^ , introduced in Refs. 14, 15 . 
The latter operators are the analogue of the l/m|k| non-Abelian-potential that was used in 
early computations where no systematic renormalization group improvement was carried out 
(see e.g. 30|). We note that the Wilson coefficient differs for quarks and for scalars 19 . 
This is, however, irrelevant for the result of the anomalous dimension and only affects 
the solution of the resulting renormalization group equations (RGB's). The diagram (e) 
contains the renormalization constant of the current which we determine in the MS scheme. 
Generically we define the renormalization constant as 

r^NLL 

Z, = 1 + ^— + ... . (36) 

e 

The resulting NLL order RGE for the current Wilson coefficients c(z/) reads 

^.|;ln[c(z.)] = 45zfL^ (37) 

Due to rotational invariance and the fact that the currents are irreducible, all compo- 
nents of any given current have the same anomalous dimension. Therefore we can compute 
the correlators for all spatial and spin indices being contracted. As an example, for the 
correlators in Fig. [2^ with the two different insertions of the potential —iV2/m'^, and after 
applying the contraction formula (fT3!l . one has to compute the three-loop diagram 

J (qi + ^) {m + ^) (q2 - qs)^ (qi + 5) 

The expression in Eq. (1551) is equivalent to the corresponding S-wave result up to the weight 
factor Nl\c\x\^\^3,\^ Piip-i qi-qs); where qi and qs are the loop-momenta on the sides of the 
diagram. The weight factor is a homogeneous polynomial of degree L of the terms q^, q| 
and qi.qa, see Eq. flT2l) . The contributions to the renormalization constant of each of the 
diagrams in Fig.l2]with the weight factor (qi.q3)'"(qi)^~(q3)^~ proper combinatorial 



17 



factors follows: 



X^NLL 



647r2 



0, 



m odd 



Fi(i,-m,|,2) + 



' 2' ■ 



vi^V)vr(^) 

, m even 



64(1 + n)7r2 
alimv) ,C 



(^(3Vi',V) + 2y<J(.)) + ^(f - c.o) 1,2) 



4 



2567r2 



(39) 



In App. |D] some useful formulae to determine these results are given. The final result for 
the anomalous dimension can then be obtained by using Eq. f|T2|) for the definition of the 
generalized Legendre polynomial. The result is remarkably simple and reads 



6z 



NLL 



IGtt 



2 c 



4 



4(2L + 1; 



647r2 



5 



LO 



(40) 



The explicit expressions for the potential Wilson coefficients for quarks and colored scalars 
can be found in Refs. [l5| and [3] (see also Ref. [21j]). For S- and P-waves the expression 
agrees with previous results, see e.g. Ref. Written in this way, the result (HOl) is valid 
for both scalar and fermion singlet currents. Differences in the running of the singlet cur- 
rents for scalars and fermions can only arise at this order from a different running of the 
Wilson coefficients for the potentials. Since at this order only the coefficient V2*'' differs for 

(s) 

quarks and for scalars, and since V2 only contributes for S-waves, the NLL order anoma- 
lous dimensions for quarks and scalars agree for L > 2. We note that the NLL anomalous 
dimension is independent of the scheme used for the currents (or the potentials) since at LL 
order the currents are not renormalized. Had we e.g. used the singlet current with cr*-'^ the 
spin-dependent potentials could in general contribute to the UV-finite terms at this order, 
but not to the NLL anomalous dimension. 



V. NLL ORDER ANOMALOUS DIMENSIONS: TRIPLET CASE 
Spin Dependent Potentials 

Before determining the NLL order anomalous dimension for the spin triplet currents, it is 
instructive to briefly discuss the form of the spin-dependent potentials from the point of view 
of scheme-dependences in the EFT in n dimensions - despite the fact that the results for 
the NLL order anomalous dimensions of the currents defined above are scheme-independent. 
From the expansion of full theory quark- ant iquark scattering diagrams without relying on 
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any a-matrix relation for n = 3 at intermediate steps, one obtains the tree-level potential 
(k = p' - p) 



'l_ p^ + p'^ 
2m2k2 



+ 



2m? 

n — 3 

Am? 



[qi \Qi Qoi] _ y y \Qi Qj^ 



2vn?\<? 



+- 



1 p2 + p'2 
k2 + 2m2k2 

n I Am?\? 



n 



1m?\? 



(41) 



where in Eqs. ( BTj) and in the following the symbol for the particle-antiparticle tensor product 
is indicated only in color space. The spin-dependent terms, displayed in the last equality, 
have been written in terms of irreducible tensors with 



p>p'\ 

-{S\ S^} 



n 



(42) 



The tree-level potential in Eq. ( 1411) is consistent with the n-dimensional expression used in 
Ref. jisl accounting for the fact that an insertion of the last two terms in the second line of 
Eq. fj^Tj) does not contribute for the current jjZL+i=i ~ 

Apart from the spin-orbit potential ~ [S**, S^]p'^p'^' which trivially reduces to ~ S.(p x p') 
for n = 3, the result does not have a very close resemblance to the well known spin-dependent 
terms in three dimensions IM 



3m2 2m2 



A(p', P) 



1 



12m'' 



T(k) 



where 



crp + (7 



ap 



A(p',p) 



. S(p^xp) 



T(k) 



3 (cTpk) {(Tapk) 
k2 



(43) 



(44) 



This is because the three-dimensional relation ® a^^ = c^a ® '^ip ~ ^^^cri 



a^p does not 



hold in this form for n ^ 3. As for the currents we are, however, free to use any scheme for 
the spin-dependent potentials as long as they are SO(n)-scalars and reduce to the known 
spin-dependent potentials for = 3. 

A suitable scheme choice can be obtained from the a structures shown in Eq. fj^Tj) . or 
from three-dimensional potentials in Eq. (1431) through the ra-dimensional replacements 



[(Tp(Tap) — 3 



-^S(p'xp) 

{(Tpk) {(Tapk) 
k2 



n 



[(Tpk) {(Tapk) 

A;2 



-2n S'^ 



T'^ik) 

k2 



(45) 



^ While any number of insertions of the spin-orbit potential ~ S^^^^T^^ ^ '1 in the correlator of two 
currents vanishes identically in n dimensions, several insertions of the tensor potential ~ 5*-'r'^(k) can 
yield non- vanishing contributions beyond NNLO. 
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We will use the scheme shown in Eq. (14 5 p for the computations of the anomalous dimen- 
sions in the following sections. This scheme is very convenient because the different po- 
tentials are total contractions of spin- and angular momentum tensors containing only one 
single irreducible representation. As a comparison, consider in contrast the scheme choice 
{(TpCTap) ~ 3(crpk)(crapk)/k^ for the tensor potential in n dimensions, which is in principle a 
viable scheme as well since it is a scalar and also reduces to the correct expression for n = 3. 
Written in terms of irreducible tensors one finds 

k^ n 

While the tensor potential based on Eq. fj^5l) is a natural minimal extension of the tensor 
potential in three dimensions, the tensor potential based on Eq. P6|) does implicitly lead 
to a change of scheme for the S^/m^ potential and the spin-independent 1/m? potential, 
which in general affects matrix elements, matching conditions and anomalous dimensions at 
non-trivial subleading order, i. e. beyond NLL order. 

We would like to point out that the spin-dependent potentials in Eqs. f H51) cannot produce 
transitions between the spin triplet currents with fully symmetric indices introduced in 
Sec. mil and other currents that transform as inequivalent irreducible representations of 
SO(n) for n 7^ 3, but that have the same quantum numbers for n = 3 (e.g. built from the 
tensors Tj"^ introduced in the Appendix. [B]). Technically this property arises because a 
given current with fully symmetric indices and quantum numbers J and L (Ej^^) always 
differs from its counterpart Ej^^ in the number of sigma matrices. Eor example, the insertion 
of a spin-orbit potential in a current with a spin structure o"*^'"*'" yields the spin tensor on 
the LHS of Eq. fl24|) . where the indices j and ^ are contracted with the initial and final 
momenta. The RHS of the same equation shows that the resulting spin structure consists 
of terms with the same number of sigma matrices as the initial current. The same holds 
trivially for the case of the potential, and also for the tensor potential, where the relevant 
relation is Eq. fl25l) . The possibility that the spin-dependent and the spin-orbit potentials 
could induce transitions between currents transforming according to different irreducible 
representations with the same J but different L is also ruled out, as the and the spin- 
orbit potentials cannot change the angular momentum L. On the other hand, the tensor 
potential can generate transitions among the Ef fully symmetric currents by two units of 
L. 



Anomalous Dimensions 



The spin-orbit and tensor potentials defined in Eq. fH51) contribute to the NLL anomalous 
dimension of the spin triplet currents through insertions in the three-loop current correlator 
diagrams in Eig.[3l A single insertion of the spin-orbit potential —iV'^^^S^'^^^ ■ T^p'p'^^ j (2m^k^) 
yields for any of the spin triplet currents defined in Eqs. (I31II32II33P 



-im 



X- 



"(E5=i(q3,o-)^-^-(za2))^ \S\S^^\ (E5=i(qi,^T)^-^-(za2)) 



?2 93 



(q? + ^) (qi - ^lY (qi + 5) (q2 - qs)^ (ql + 5) 



(47) 
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FIG. 3: Three-loop current correlator diagrams with insertions of the spin-dependent potentials 
for the computation of the NLL anomalous dimension of the currents. Combinatorial factors are 
suppressed. 



where the proper contraction of spinor indices is in analogy to Eqs. ( l30|MIl . Note that the 
correlator of two spin triplet currents with the same J but different L vanishes, as the spin- 
orbit potential commutes with the squared angular momentum operator L^. To compute 
the contraction in the numerator of Eq. P7|) it is useful to perform the shift (ITSl) in the 
qi integration. This reduces all non-trivial contractions in the numerator to the form of 
Eq. (fT5|) after use of the relations (fTOj) . The computation then reduces to 



-im 



with 



|qi|^|q2| IqsI^^^^lI^, qi-q2) (i^L+il?^, q2-q3) - q2- 

D-c^, D-cis (qf + 6) (q^ - q^)^ (qi + 6) (q^ - qs)^ (qi + S) 



-L{L + l){2L + n)^^, J = L + 1 
cj,L=\ L{n-2){L + n-2) , J = L 

L + n-2, J = L-l 



qs. 



(4J 



(49) 



The calculation of the three-loop current correlator with a tensor potential insertion 

can be carried out similarly. For the case where the triplet 
currents at the sides have the same L the correlator reads 

-im N, Vi') Vt\-2n) j D"qi D"q2 ^"qs 

(rfHq3,'^)^"-"^(^^2))^ 5*^T*^(q3-q2) (Ffi(qi,fT)^i-'^(za2)) 



Tr 



(q? + 5) (qi - q2)2 (qi + 5) (q2 - qa)^ (qi + 5) 



-imN,V'f^vi:'^NLCj,L'Ti\l] / D'^qiL>"q2L>"q3 



|qil^|q3l^^L(^, qi-q2) xl(^; q2, qs) 

(q? + 5) (qi - q2)2 (qi + 5) (qi + 5) 



(50) 



with 



XL{n] p,qj 



[n 



— n ■ 



2)(L + n-2)Pi(n,p.q) 

{n-2){L + n-2) |p||q| 
2L + n - 2 (p - q)2 



PL-i{n, p.q) - PL+i{n, p.q) 



(51) 
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FIG. 4: Three-loop current correlator diagram between the currents with orbital angular momen- 
tum L and L' = L ±2 and a tensor potential insertion. 



and 




2L+n-4 



(52) 



The divergent part of Eqs. (H8l) and ( ]50l) can be extracted from the integrals listed in 
Appendix D. The insertion of a potential is identical to that of the V2/m^ potential, 
Eq. (|38|) . and the result can be read off from Eq. fBOl) . The result for the renormalization 
constant of the currents from the diagrams in Fig. [3] including combinatorial factors thus 
reads {S = 1) 



c,spin 



.NLL _ V^^^(^)V1^^(^) J(J+l)-L(L + l)-2 ,^ 

- 256^^^ (2L + 1)L(L + 1) " 

TTT^ dLoS{S + l), (53) 



647r2 2L + 1 ' 647r2 



1 



J(2J+1) ' J — 1/ + 1 

^,^=<l Tf;^' J = L (54) 

J = L-1 



(J+1)(2J+1) 

Explicit expressions for the Wilson coefficients of the spin-dependent potentials can be found 



in Ref. [l6|. For "^^^^Lj = ^Si the result agrees with the known results from Refs. ll|, [l5|, 
18l . |22j, and for '^^~^^Lj = ^Pi q the renormalization constant at the matching scale z/ = 1 is 



consistent with the IR-divergences in the hard region obtained from the threshold expansion 



in Ref. 47 



The tensor potential can also change the orbital angular momentum L by two units while 
conserving the total angular momentum J. ^ Thus transitions between spin triplet currents 
with orbital angular momenta L = J ^ 1 and L' = J±1 = L±2 driven by diagrams such 
as Fig. m with a tensor potential insertion do not vanish and could induce a mixing between 
the currents through renormalization. The expression for the three-loop correlator of Fig. HI 
where the tensor potential is inserted next to the annihilation current on the RHS, reads 
now 



^ Note, however, that the spin of the fermion-antifermion system (S=0,1) cannot be changed by the S = 2 
tensor S^^ . Changes of the angular momentum L by one unit are not aUowed due to parity. 
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NM'^Vi'\2nJNj+,)Ti[l] 

n" Iqil^^lqal-^'^il^, qi-q2) q2, qa) 

Z; qiZ; qzD qs 



(qf + 5) (qi - q2)2 (ql + S) (q^ - qs)^ (qi + 6) 



(55) 
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XJ,L[n; P,qj 



(56) 

The divergent contributions of Eq. fl55|) can be obtained from the results for the relevant 
integrals given in Appendix D. We find that the UV divergences from the sum of both 
insertions of the tensor potential cancel and that there is no mixing between the and 
the jj=L+i currents at NLL order. 



VI. SOLUTIONS AND DISCUSSIONS 
Scalar Production Currents 

For scalar fields, only the spin-independent contributions shown in Eq. (I40l) have to be 
taken into account and we can write the RGE for the Wilson coefficient of the current with 
quantum number L in terms of the RGE's of the S- and P-wave current Wilson coefficients 
which were determined in Ref. 



d 



d\n.v 



NLL 

c 



'^MTr^ln[c5(z^)] + 
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The solution is thus found immediately and reads 



'LO) 



dlnu 



ln[cp(z/)] 



In 
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+ i?4 Ttasim) 



1 _ ;21-13Ca/{6/3o) 



Bq -nasim) 
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where z = as{mi')/as{m), w = asimv"^) / as{mv) , and 



(57) 
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Bi = Slo h + 



2L 



'I - 6lo) di , 



(59) 
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and with 1^ 



dn 



Cf[CaCf{9Ca - WOCf) - f3o{2QCl + 19CfCa - 32C^) 

Cj, [CAjOCA - IOOCf) - 6Po{3Ca - 4C^) ] 
2/32(6/5o - 13Ca) 

24CUnCA-3f3o){5CA + 8CF) 
13CA{Qf3o - ISCa^ 

8CFiCA + CF)iCA + 2CF) 

Cf{Ca + 2Cf) 



3/?o 

8CACF{CA + iCF) 

9(31 



(60) 



In Fig. Owe have displayed the NLL running of the scalar currents exemplarily for L = 0, 1 
and 4 normalized to their matching values. For the input parameters we have chosen two 
different values for the heavy scalar mass, m = 220 GeV (solid lines) m = 500 GeV (dashed 
lines), and as{mz) = 0.118, taking leading-logarithmic running for with Uf = 5 active 
massless quark flavors and no active massless squarks (n^ = 0). As is already visible from 
the form of Eq. (1571) . we find that the i/- variation of the currents becomes weaker for larger 
angular momenta. While the S-wave coefficient increases by about 6% for u ^ as compared 
to the matching value ai u = 1, the P-wave coefficient only increases by less than 3%. For 
L < 4 the maximal relative variation is already below 1%. It is also conspicuous that for 
any L the maximum slightly decreases with the heavy scalar mass and also moves towards 
smaller values of u. As already discussed in Ref. [l9|, the latter feature can be understood 
qualitatively from the fact that the average velocity (v) ~ as{mas) decreases with the heavy 
scalar mass m. 



Fermion Production Currents 

The RGE of Eq. fl57|) for the Wilson coefficients of the currents which produce a pair of 
quarks with quantum numbers J, L and S is obtained at NLL by adding the anomalous di- 
mensions from the diagrams with spin independent and spin dependent potentials (Eqs. fHDj) 
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and (|53|) respectively): 
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(61) 



We have singled out the known RGE for the S-wave currents UJ, llJ, l22| , Cjq, that is relevant 
for 7* — >■ qq production in the spin triplet case, and for 77 — qq in the spin singlet case. 
Solving Eq. (16T|) yields {z = as{mh')/as{m),w = asimu"^) / asimi')) 



In 
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(62) 
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The coefficients Oj for the running of the S-wave Wilson coefficients Cjq ' read [U, |22 
02 = 62 , 04 = 64 , ao = &o , 



(63) 
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(64) 



where = 2 in the triplet case and = for the singlet case. 

In Fig. and b we have displayed the NLL running of the quark currents for spin singlets 
and triplets, respectively, for a number of cases. For the heavy quark mass we have chosen 
m = 175 GeV, and for the strong coupling we employed the convention used in the discussion 
of the scalar currents. For the singlet cases we have shown the evolution for L = 0, 1,2,4. 
The dependence of the evolution of the currents on the value of L is very similar to the 
scalar currents. For the triplet cases we have displayed the evolution for the cases L = 0, 1, 3 
and all possible values of J. Here the dependence of evolution on the choice of the angular 
momentum L (for a fixed J value) is similar to the spin singlet cases. But we also find that 
the evolution for a fixed value of L is stronger for the smaller J values. 



VII. COMMENTS ON SCHEME DEPENDENT RESULTS 



In this section we comment on scheme dependences of a number of results for nonrel- 
ativistic currents that can be found in the literature. The current that has been stud- 
ied most extensively is the vector current 'ip'^^CT(ia-))'}(*_^, which governs heavy quark pair 
production in e'^e 
e+e" ^ UH E 



annihilation [3l|, |50|, l51| or the large Higgs energy endpoint region in 
33| . At NNLL order the matching condition and the anomalous dimen- 
sion depend on the scheme used for the spin-dependent interactions and potentials. Using 



dimensional regularization such computations were carried out in Refs. (I4j . Il5l . |29| . |4J, |45 
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V V 

FIG. 6: Left panel: NLL running of the normalized Wilson coefficient for quarks in a spin singlet 
configuration for L = 0, 1, 2 and 4 for m = 175 GeV. Right panel: NLL running for quarks in spin 
triplet configuration for L = 0, 1,3 and all possible J values for m = 175 GeV. 



In Refs. 14, 15, 2^, the NNLL order matching condition of the vector current 
ipp(T{ia2)x*-r, was determined using three-dimensional a structures for the spin- dependent 



potentials, as shown in Eq. fH3l) . and three-dimensional a algebra. Although the use of the 
three-dimensional a algebra is a priori inconsistent, the results for the NNLL order matching 
conditions derived in Refs. [3, 15, 2^ nevertheless represent a viable scheme. This scheme 
is related to the spin- dependent potentials in Eqs. (14T!) or (HSl), with strict n-dimensional 



a algebra, through multiplicative n-dependent factors. For the vector current 
this scheme has, for the tree level potentials, the form 



1 , P^ + P^^ 
k2 2m2k2 
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—A--1 
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k2 



12 



P" + P 



+ (l + e) 



3 m? 



(65) 



and can be derived from the entries in Tab. [TTl We have dropped the spin-orbit and tensor 
potentials that were displayed in Eqs. (HTl) or (145!) because they do not contribute for the 
vector current '?/'^<T(i(T2)xlp. The first relation in (1651) allows to compute the difference 
between the NNLL matching condition obtained in Ref. [15] and the scheme where the NNLL 
matching condition is identified with the contributions from the hard region in the threshold 



expansion 



4J, 45], which implies the use of the scheme for the potentials shown in Eq. (14 ip 
[see also Ref. 52|)- These two schemes only differ in the treatment of the spin- dependent 



potentials. Using the contribution to the NLL anomalous dimension from the spin- dependent 
potential Vs at the hard matching scale from Eq. (|53 



m 



)/6 0, [17 



^ NLL 



-Vi^^(l)Vr^(l)/327r2 



one finds that the relation between the NNLL matching condition 

'1), and the one based on 



based on the three-dimensional a algebra from Ref. 15|], c 
the threshold expansion scheme, c^nll'v (-'-)' reads 



„3dim 
-NNLL,' 
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'-NNLL,vl-'-y 



thr exp 

~ %NLL,v 



_ Aa^NLL _ threxp /.n _ ^ fi2 2 
-|^2 c.spin — %NLL,vl-'-J ^^^F^^s^ 



ml 



(66) 



which agrees with the results given in Refs. |l5l . |4J, |45j. It is also straightforward to com- 
pute such a scheme relation for the spin-dependent non-mixing contributions of the NNLL 
anomalous dimension of the vector current that was also computed in Ref. 151 . 
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Scheme relations for the NNLL matching conditions analogous to Eq. (166|) . connecting 
the scheme from Ref. [l5| with the one based on the threshold expansion can also be derived 



for the ^Pi axial- vector ( J*^ = 
pseudo-scalar {Jpg = ilj^{ia2)x^-. 



„) currents: 



^Po scalar (J, = ^l^a^{ia2)X-p), and ^Sq 



^3 dim 
-'NNLL,av 
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^3 dim 
-NNLL 



.s(l) 



4niZs(1) -^C^pOilirn) 



„3 dim 
-NNLL.ps 



(1) 



thrcxp /-| \ 
^NLL.psUJ 



(67) 

(68) 
(69) 



The respective NNLL matching conditions in the scheme based on the threshold expansion 
were computed in Ref. [47| . 

Finally, we would like to make a few comments on the conventions for 75. In previ- 
ous NRQCD literature based on dimensional regularization, only the totally anticommuting 



47 , which leads to the nonrela- 



version with the form 75 = (jq) has been considered 

tivistic currents V'p(^'^2)Xp and ip^j^a^\ia2)Xp ^^i the pseudo-scalar ^ip^ip and axial- vector 
"07*75 = 1, . . . ,n) full theory currents, respectively. 



In particular, in Ref. [47l] the 
NNLL (two loop) matching conditions for these pseudo-scalar and axial-vector currents 
were identified with the hard contributions obtained in the threshold expansion using the 
form 75 = Z7''7^7^7^ = — «7o7i7273 for triangle graphs and the totally anticommuting 75 
for the other diagrams. Such an approach is viable as long as the hard contributions from 
triangle graphs only lead to IR-finite terms. However, from the point of view of the full 
theory it can be important or even mandatory to employ only the fully consistent defi- 
nition 75 = Z7°7^7^7^. The latter form for 75 also leads to a different from of the EFT 
currents from the nonrelativistic expansion. Writing 75 = ^e^^^'j^Y'j^'-f^ 



and 

where e^^^ is the usual three-dimensional e-tensor if all indices assume 
values 1,2 or 3, and identically zero otherwise, one obtains for the leading order expansion 



7 75 



^7*75^ 



m 



{^cr2)X-p 



(70) 
(71) 



The current in Eq. (!70!) was discussed in Sec. Illll and the one in Eq. (17T!) is a simple example 
for the alternative ^L^ currents involving the tensor F^"^ shown in Eq. flC3p . The EFT 
currents on the RHS's of Eqs. fl7U]l and (17X1) represent a viable scheme and are equivalent 
to the currents above for n = 3, but inequivalent otherwise. So at NNLL order their 
respective matching conditions and anomalous dimensions differ. Using 75 = Z7°7^7^7^ in 
the full theory, the 2-loop hard contributions from the threshold expansion give the NNLL 
matching conditions of the currents in Eqs. (ITOll and (17T|) . 



VIII. CONCLUSION 

The construction of non-relativistic interpolating currents for describing the production 
and annihilation of a color singlet heavy particle-antiparticle pair with arbitrary quantum 
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numbers ^^'^^Lj in n = 3 — 2e dimensions has been addressed in this work. Arbitrary 
angular momentum configurations are accounted for by the generahzation of the spherical 
harmonics in n dimensions, for which we have provided a representation in terms of cartesian 
coordinates. The use of the n-dimensional spherical harmonics is required to maintain SO(n) 
rotational invariance in calculations within dimensional rcgularization. 

Similarly, a consistent description of the spin of a fermion-antifermion pair requires the 
treatment of Pauli cr-matrices in n dimensions. The antisymmetrized product of cr-matrices 
are the proper irreducible tensors with respect to SO(n) to build interpolating currents 
describing an arbitrary angular-spin state of a fermion-antifermion pair. In n dimensions 
the form of irreducible currents in SO(n) which reduce to a specific angular-spin state in 
n = 3 is not unique. This is related to the existence of evanescent operators that can 
be built from the nonrelativistic field operators and the cr-matrices in n dimensions. Such 
evanescent structures occur even in simple standard processes, and in general need to be 
taken into account consistently at subleading order for the renormalization process. The 
specific choice of a basis for the physical currents defines a specific renormalization scheme. 
Similar scheme choices are possible for the spin-dependent interactions and potentials. In 
this work we have discussed different versions for spin-singlet and spin-triplet currents, and 
also for the spin-dependent potential, which are inequivalent in general n dimensions, but 
equivalent for n = 3. 

We have discussed the latter issues in the framework of NRQCD and have shown that 
at NNLL order transitions between currents containing different irreducible a structures 
cannot occur. Since for a specific angular-spin state the corresponding leading order currents 
do not have an anomalous dimension at leading-logarithmic order, the scheme- dependence 
introduced by specific choices for the currents and spin- dependent potentials in NRQCD 
only affects NNLL matching conditions and anomalous dimensions. The relation between 
NNLL matching conditions of a number of currents for schemes used in the literature have 
been determined. 

Finally, we have determined the NLL order anomalous dimensions for the leading in- 
terpolating currents for quarks-antiquark and colored scalar-antiscalar pairs for arbitrary 
spin and angular momentum (J, L, S) quantum numbers. We have shown by explicit com- 
putation that the spin-dependent potential do not cause any mixing between currents with 
different L at NLL order. An important property of the solution of the respective anomalous 
dimensions is that the scale variation of the currents becomes in general weaker for higher 
angular momentum L. 
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APPENDIX A: FERMION PAIR PRODUCTION FROM 77 COLLISIONS 



The 77 — > // amplitude at lowest order in a is given by the graph in Fig. [7] plus the 
one with permutted photons. Let us write the production amplitude up to second order 
in the three-momentum of the fermions in the cm. frame p. The outgoing momenta are 
p = [E^ p), and p' = (i?, — p) with E = a/ + p^ = m + ^ + . . . , and the three- momenta 
of the photons are ki = — k2 = k. We shall use Pauli spinors ip and x fo^' the fermion and 
the antifermion respectively, 



E + m 
2E 



m+E V^P 



^ + ^ £^ i^^2)X 



2E 



m+E \""^J^p 



(Al) 



and pick a gauge where the photon polarizations are purely transverse, €{■ ki = 0, e° = 0. 
The result for the full amplitude in Fig. [7] including the crossed graph is 

1 



— — m 



i2v{-v) + {I ^ 2} 



(A2) 



and expanding in p gives 



iM 



m 



ei ■ €2)^1 7 



m(eie2) (m + 7^') + M ^e!, 7*^' 



m 



-p) + .... 



(A3) 



The spinor structures in Eq. flA3p can be further expanded in p by using the leading or- 
der matching relations between full theory currents and the non-relativistic ones shown in 
Table E 



where 



+ (p' + pi a' - ^5''Vp)) {^(y2)X-p + C>(pV^' 



(A4) 



-.2 



/ J J 



m 



/ 



2m? 



€162) 



+ e\ 4 + 4^2 



Thus in fermion production from 77 collisions the currents j|*£^Q (^5'o), jf_i (^-Po) and (jf+i)*-' 
(^P2) are involved. 



(A5) 



APPENDIX B: FERMION- ANTIFERMION PAIR ANNIHILATION INTO 87 

The decay amplitude // — > 37 to lowest order in a is given by the graph in Fig.[8]summed 
over the 3! permutations of the photons. The amplitude of the graph shown is 

,M = -.e^Q?.(-p)fe _^,/ /.^3|^A^P). (Bl) 
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k2 ^AAA> 



pi 



FIG. 7: Fermion pair production from two photons. 



-\AAAy /ci 
'^AAA/ h 
-\AAAy A:3 



FIG. 8: Three photon annihilation graph. The 3! permutations of the photons have to be consid- 
ered. 

We want to compute the leading term in the expansion of the three-momentum of the 
fermions in the cm. frame p. We can then work at threshold and the spinors reduce to 



«(p) 







(B2) 





The leading term of the amplitude corresponding to the photon ordering in Fig. [8] then reads 



Am? 



(cre3)(crfc3)(cre2)(crA;i)(crei) + (cre3)(cre2)(crei) 



p ' 



(B3) 



where fc, = ki/uoi are the normalized three-momenta of the photons. If we add the graphs 
with different ordering of the photons and write the products of Pauli matrices in terms of 
the antisymmetric sigma tensors we obtain for the // — > 87 amplitude: 



iMo = i 



-,3/n3 

2 X-p {-^<^2) [ —(-I Kl £2 ^3 ^ 



m 



^3 



/(I, 2, 3) e{ + g{l, 2, 3) ki + {1 ^ 2} + {1 ^ 3} , (B4) 
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FIG. 9: Resolution of the identity in terms of Young operators of the symmetric group Sl+i 
accounting also for the proper Hook factors [HsI]. We have only shown the Young patterns which 
give a nonzero result when acting on the indices of J'*^ "*^c7*^+i. We shall use [L + 1] to denote 
the first pattern, and [L, 1] for the pattern with L boxes in the first row and one in the second. 
There are L standard tableaux for the pattern [i^, 1]. We adopt the convention that symmetrization 
horizontally top-to-down is followed by antisymmetrization vertically left-to-right. 

where 

/(I, 2, 3) = (l + (^2^3) - (fcifca) - (fcifc2)) (€263) 

+ (fcie3) (^362) + (fcie2) ik2es) - (^62) {k2^z) 

g{l, 2, 3) = ((fcsei) + (^261)) (€263) - {e^e^) (^362) - (£162) (^263) • (B5) 

Note that the a*^*^*^ terms cancel when summing over permutations of the photons. 

APPENDIX C: TENSOR DECOMPOSITIONS 

1. X^i-iL (jiL+l 

The space of traceless tensors of rank L + 1 formed from the product of the T*^ ' *^(p) 
and 0"*-^+^, ^ 2"«i---«icr«i+i^^ is obtained by removing the traces of 



fc=i 

L+l 



2 

(L + n-3)(2L + n-2)^^ ^l-iIp,'^) , (^1) 



with the tensors Tf^} defined in Eq. (1321) . This space has dimension ul ■ n — ni-i and 
is not irreducible. From the resolution of the identity in terms of Young operators of the 
symmetric group Sl+i [4o'] (see Fig. [9]), which act on the indices of the tensors, we can 
decompose ^ T*^'' *^cr*^+i) into traceless tensors of a given symmetry type: 

k=2 
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The first term in Eq. (lC2p results from tlie application of the first tableau in Fig. [9] (the 
pattern [L + 1]) to the tensor T*^'' '-^cr*^+i), which yields a totally symmetric tensor. This 
is the tensor used to define the J = L + 1 triplet current in Sec. IIIIl and its explicit form is 
given in Eq. fl3ip . The L terms in the second an third lines of Eq. (lC2p correspond to the L 
Young tableaux with the pattern [L, 1] shown in Fig. O These tensors read 

. L+l 

= T*i-^'"-*^+i a^'" V (5*'=*™ T*i-^'=-^™-^^+iV^ - Ui ^ , (C3) 

T , -\- n — ' ^ 



L + n-3 



and have dimension (L > 1) 

V{n + L-l) L{2L + n-2) 

d[L,i] = riL-n- nL_i - n^+i = — — — — , C4 

1 (n — 2)1 (L + 2) L + n — 6 

which for n = 3 reduces to 2L + 1. Indeed, the tensors corresponding to the pattern [L] in 
Eq. ( l33l) and those corresponding to the pattern [L, 1] in Eq. (EH]), although transforming 
according to different irreducible representations of SO(n), become equivalent in = 3, and 
thus represent both appropriate candidates for defining the triplet current with quantum 
numbers J = L in dimensional regularization. This justifies the quantum numbers used to 
label the currents in Eq. (]C3p . 

As the tensor (( T*^'"*^cr*^+i) is already symmetric on the ii. . Al indices, the subspaces 
which result from the action of the different [L, 1] tableaux are all equivalent. Therefore the 
space ( T*i'"*-^(T*-^+i) decomposes into a subspace of totally symmetric tensors, of dimension 
rii+i, and a subspace of tensors which transform equivalently to the irreducible representa- 
tion of SO(n) defined by the symmetry pattern [L, 1], of dimension d\^L,i]- Thus the second 
and third lines in Eq. (lC2p should not be understood as a tensor decomposition into orthog- 
onal subspaces, but as a convenient way to write the components of a tensor of this subspace 
in terms of tensors of the symmetry type given by the Young tableaux. 

The results in Eqs. fICip and flC2p show how the reducible tensor = 
2iii...iz,^p^^«£+i^^^^^^^^ pg^j^ fully written in terms of the currents defined from the 

tensors Ff , Ff"^ and Ff^^ The latter current already emerges in Eq. (ICip when elimi- 
nating the traces from T'i'"*^o"*^+i. 
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FIG. 10: Resolution of the identity in terms of Young operators of the symmetric group S'l+2- We 
have only shown the Young patterns which give a nonzero result when acting on the indices of 
2iji...jlu«l+i«l+2_ There are in total 2(L — 1) tableaux in the second line. 

2_ J^il-.-lL fjiL+liL+2 



The reduction of the tensor T*i "*^(p)cr*^+i*^+2 can be carried out in a similar way. The 



traceless tensor (| r*^ "*^(p)cr*^+i*^+2 ), of dimension Ul 
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ni_i ■ n + 711-2, reads 
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The latter is a tensor which transforms according to the representation defined by the pattern 
[L, 1], so for n = 3 it is equivalent to the totally symmetric tensor r|^}(p, ctY^-^l-i defined 
in Eq. fl32|) with L — 1 indices. 

The space ^ T*^ ' *^(p)(T*^+^*-^+2 ) can be written in terms of tensors with the symmetry 
type given by the standard Young tableaux which appear in the resolution of the identity 
of the symmetric group 5*^+2 shown in Fig. [TD) 



'rph...iL^iL+liL+2\ = rf=|(p, Cr)!^-."^+|*^+^ — FflKp, (t)!'-."^+/ + . . . (C7) 

/ L -\- 2 -'^+-^^ ^[mL+i] L + 2 'lmL+2] ^ ' 
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The first and second terms in Eq. (1C7P correspond to the action of the two Young tableaux 
with the pattern [L + 1, 1] in Fig. [TOl The explicit form of these tensors read 



L+l 

"r5'=l/' y2---2L + l _ rpi2---iL+l _|_ \ ^ rp il-.^k ■■ -iL+l ^ikiL+2 



k=2 

L+l 



k<i 
fc=2 



L+l 

k,j=2 



L+l 



- ifl + /2) Yl T*^-^'=-*^+lV"^+2 - {^1 ^ tL+2} 



k=2 

f = ^ f = 2L + n 

■'^ 2L + n-2 ' ■'^ {L + n-2){2L + n-2) ' ^ ' 

The dots in Eq. (1C7P stand for the tensors associated to the rest of the Young tableaux in 
Fig.[TOl We do not give their explicit form because they are evanescent and vanish for n = 3, 
as the traceless tensors corresponding to Young patterns in which the sum of the lengths 
of the first two columns exceeds n are identically zero ^. Eqs. f ICSP and flC7l) give the 
decomposition of the tensor _B*i" *i+2 _ T*^'' *^(p)(T*^+^*^+2(i(j2)xlp into the currents built 

from the irreducible tensors rf^{ and Ff^^, that have fixed quantum numbers in n = 3. 

The three new spin triplet currents ipp^j^^ {'i<^2)X-p introduced above, Eqs. ( lC3fC6IIC8p . 
have 2J + 1 independent components in n = 3 dimensions and differ from the basis of 
currents defined in Sec. IIIII in the symmetry patterns and the number of indices. The 
following relations among the tensors Fj^^ and Fj"^ with the same quantum numbers hold 
for n = 3: 

ff=}(p,^)i;;^^| = ^^6-^^^-^Ff=l(p,ar-^-S 



For the case L = 1, the third diagram in Fig. [10] does not vanish in n = 3, but yields a totally antisym- 
metric tensor. 



35 



For completeness, we also give the total contraction of the Tj ^ tensors: 



Tr 



L+l ' 



p,cr 



,«2 ^L+l -p5=l 



«2 

iliL+2] 



Tr 
Tr 



22-. .^L F'5=l/' 



2(n - 2)(L + 2)(2L + + n - 1) 
2L + n - 2 

xTi,„(p,q)Tr[l], 
, L + n-2 



2(n 



v«2---«L-l -p5=l 



p5=lt/ 



' L + n-3 
2(n - 2) Ti,„(p, q)Tr [1 



7'L,n(p,q)Tr [1] , 



(CIO) 



To conclude this section we mention that the guidelines of constructing nonrelativistic cur- 
rents that are irreducible under SO(n), as shown above, can be generalized in a straightfor- 
ward way to currents containing o"*!- -*'" with m > 2. For larger m, the number of evenescent 
currents will increase. A very helpful rule to identify these currents is the rule that traceless 
tensors corresponding to Young patterns in which the sum of the lengths of the first two 
columns exceeds n are identically zero (4o| . 



APPENDIX D: USEFUL INTEGRALS 



We list in this appendix the results for the relevant integrals needed for the determination 
of the UV-divergences of the three-loop correlators in Secs. lIVI andlVl Note that Greek letters 
a, /3, 7, . . . refer to real numbers while k, r denote positive integers. 

We define the following integrals: 



I^'^[«'A7] 



(p-q)' 



(2^ (2^ [p2 + 5]" [(p - q)2]/3 [q2 + 5]' 



(Dl) 



l(^')[a,/3,7,P,a] = 



rf"p d'^i rf"q (p.q)*^ 

(27r)" (27r)- (27r)" [p2 + 5]" [(p - £)Y + S]^ [{i - q^Y [q^ + 5]- ' 

(D2) 



J[a,/?,7,p,(T] 



rf"p d"q 



1 



(27r)" (27r)" (27r)" [p2 + 5]° [(p - £)2]/3 [£2p [(£ _ q)2]p [q2 + ^ja 
For the spin-independent potentials the required UV-divergent terms are: 

I^^^[l'^'lW^(-l)-.^(^'-4'2)^ + O(e0), 



■ (D3) 



(D4) 



I(^)[1,2-|,1K 



-1)'=+SFi(^,-A;,-,2 



3 
2' 



+ 0(e°) 



(D5) 



l(')[l,l,0,l,l]div 



(D6) 
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1 r^2A:+i 

,(-.|l,0.1.1.1K = -^^ + O(eV 
The spin- dependent contributions require: 

J|l.l.-*.l.lW = (-l)'«^^ + 0(.»), (D7) 

J[l, -k, r, 1, IW - i-^r \^,,\l4ri2r) 256^ + ^^^^^ ' ^^^^ 



i [1,1, /c,l,ljd„ - ^^^^ 3) 256vr3e+^^'^ 



(D9) 



IW[1,-A;,1,1,1W = (-I)'^'t-T V + (DIO) 

1 + fc TT-^e 
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